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A brief review is given of attempts to understand tlie energies of four-quark systems 
calculated on a lattice in terms of nuclear-physics-inspired many-body techniques 
^ ' involving interquark potentials. Results are given for the next stage of this study 

1/^ , where the wavefunctions of heavy-light mesons are also calculated on a lattice. 

o, 

On , Over the past few years the authors have been measuring on a lattice the 

3^ ' energi£a.|Df various four-quark systems. In the original papers (see for example 

^— V , Refs. aaB) the four quarks involved were all considered to be infinitely heavy. 

The resultant energies could then be reasonably well understood in terms of a 
C^ , many-body nuclear-physics-inspired approach involving interquark potentials 

- provided there was introduced a four-quark term similar to a form factor. 
Oh, Neglecting such a factor consistently led to an overestimate of the binding. 

Later in Ref. □ a method was developed for treating on a lattice two quark 
systems, where one of the quarks was a light quark i.e. the case of heavy-light 
mesons Qq. In that paper the authors concentrated on measuring the S-, P-, 
D- and F-wave energies. 



H I Returning to the four-quark system, in Refs. U the energies of the Q^q^ 

system were calculated using the same techniques that proved successful in 
Ref. U for the ha^ic Qq case. In addition to the presence of light quarks, the 
works of Refs. uEI had two other improvements compared with Refs. uHu: 
i) The gauge group used was SU(3) and not SU(2). 
ii) The lattice configurations were unquenched. 

In Ref. u the earlier nuclear-physics-inspired approach in terms of in- 
terquark potentials was.^extended to the Q^q'^ case. This required fitting first 
the Qq energies of Ref. □ to extract an effective light-quark mass of about 400 
MeV. The main conclusion from this work was that the Q^(f data could not 
be understood in temw of purely two-quark potentials and, as in the earlier 
static case of Refs. EJEIu, a four-quark form factor was necessary. 



Most of the above work has been devoted to the energies of the various 
quark systems - the exception being Ref. Q where flux-tube structures were 
measured. Now we are working on a lattice measurement of the radial wave- 
functions of a single heavy-light meson. These wavefunctions consist of the 
distribution of the light quark and the colour field components around the 
static quark. The light quark wavefunctions of the ground state and some ex- 
cited states are being measured. Such wavefunctions have not been measured 
before and are of relevance to various phenomenological attempts to reproduce 
meson-decays and scattering of mesons. These include e.g. bag models and 
semirelativistic Schrodinger and Blankenbecler-Sugar equations. 

The actual wavefunction measurement is based on the light-quark propa- 
gators Gij of Ref. □. For a measurement of the Qq energies only one Gij enters 
in the 2-point correlation as, essentially, C2(i) = Xlii ^ji^ij where Uij is the 
static quark propagator represented by a straight line of gauge links from point 
i to point J in a different time slice. However, for the light quark wavefunc- 
tion measurement two such operators arise giving a 3-point correlation of the 
form C^it^r) = J2iji GjiOGiiUij, where the site I is constrained to be within 
r spacings from the i,j space coordinates. Here we use the local operators 
O = 74 and 1, which are probing respectively the light quark charge and mat- 
ter distributions at a distance r from the heavy quark. The latter are defined 
as {G3{t,r)/C2{t)). The result of fitting these distributions by F^{0), where 
F — Aexp(— r/rg), is given in Table ll|. There it is seen that the charge dis- 
tribution has a considerably longer range than that of the matter. Summing 
over the charge distribution should give the charge of the quark. With the 
present normalisation this should be wl on a lattice and, within the expected 
accuracy, this is indeed the case, when the sumJs carried out directly on the 
lattice - see the column DSum. As discussed inH the sumrulc for — 1 has a 
less direct interpretation. 

The authors wish to thank the Center for Scientific Computing in Finland 
for their cooperation in making these studies possible. 



Table 1: Parameters for fitting the charge (74) and matter (1) distributions with F'^ , where 
F = Aexp(— r/ro). Dsum refers to a direct lattice estimate of the sum of F^ . 



Operator(O) 


ra/a 


A 


DSum 


74 (Charge) 
1 (Matter) 


1.56(2) 
1.15(5) 


0.45(1) 
0.46(2) 


1.12(5) 
0.25(5) 
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Figure 1: The correlation (Cs{t,r)/C2{t)) as a function of r in lattice units: a) Charge and 
b) Matter. Solid(dotted) for t = 8(10). 



